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We have investigated the spin fluctuations in the langasite compound Ba3NbFe3Si20i4 in both the 
ordered state and as a function of temperature. The low temperature magnetic structure is defined 
by a spiral phase characterized by magnetic Bragg peaks at g=(0,0,r ~ 1/7) onset at Tjv=27 K as 
previously reported in Ref. [2l]. The nature of the fiuctuations and temperature dependence of the 
order parameter is consistent with a classical second order phase transition for a two dimensional 
triangular antiferromagnet. We will show that the physical properties and energy scales including the 
ordering wavevector, Curie- Weiss temperature, and the spin-waves can be explained through the use 
of only symmetric exchange constants without the need for the Dzyaloshinskii-Moriya interaction. 
This is accomplished through a set of "helical" exchange pathways along the c direction imposed by 
the chiral crystal structure and naturally explains the magnetic diffuse scattering which displays a 
strong vector chirality up to high temperatures well above the ordering temperature. This illustrates 
a strong coupling between magnetic and crystalline chirality in this compound. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 



I. INTRODUCTION 

In geometrically frustrated magnets, each spin cannot 
satisfy all pair wise interactions as a result of the crystal 
symmetry and therefore remain disordered to temper- 
atures well below the Curie- Weiss temperature (Qcw) 
where magnetic order is expected. ^"■^ The simplest geo- 
metrically frustrated system is arguably the two dimen- 
sional triangular lattice where all interactions cannot be 
satisfied resulting in large degeneracies.^ 

Low-spin (or quantum) triangular antiferromagnets 
with 5 < 1 do not show long-range magnetic order 
because of the combination of strong quantum fluctua- 
tions and geometrical frustration and a number of in- 
teresting phases have been observed and proposed in- 
cluding resonating valence bond states and spin liquids.^ 
Examples of two-dimensional (2D) triangular antifer- 
romagnets include the CS2CUCI4 (S=l/2), k-(BEDT- 
TTF)2Cu2(CN)3 (S=l/2), ZnCu3(OH)6Cl2 (Kagome) 
and NiGa2S4 (8=1).^^"-^^ All of these systems display 
anomalous magnetic properties well below Qcw and il- 
lustrate the dramatic effects that crystal symmetry im- 
posed degeneracy, or geometric frustration, can have on 
the magnetic structures and excitations. 



Large S systems are interesting to investigate as quan- 
tum fluctuations are completely suppressed and therefore 
reveal only the classical effects of frustration. An inter- 
esting example is the case of the anisotropic triangular 
lattice Q:-NaMn02 where the S — 2 spins, with domi- 
nant nearest neighbor coupling, relieve the effects of frus- 
tration by reducing the dimensionality and behave like 
one-dimensional chainsi^i More complex systems include 
(Tb,Ho)Mn03 (Refs. MM) and RbFe(Mo04) (Ref. M) 
where frustration results in spiral phases and coupling 
between magnetic and ferroelectric order parameters.— 
For some of these systems, the loss of inversion symme- 
try in the ordered magnetic state allows antisymmetric 
exchange (Dzyaloshinskii-Moriya interaction) which has 
been proposed to be very important to multiferroicity. 
R is therefore important to investigate the property of 
chirality and its relationship to the crystal structure es- 
pecially in systems that are non-centrosymmetric already 
in the paramagnetic statci^ 

An interesting class of compounds is the langasite se- 
ries of materials which consist of isolated triangular lat- 
tices (or trimers) on a hexagonal lattice. Attention to this 
system started with La3Ga5SiOi4 (LGS) which displays 
favorable piezoelectric properties and strong electrome- 



chanical coupling constants. ^^ LGS has a Ca3Ga2Ge40i4 
type structure with space group P321 and this structure 
type can accomodate a number of different cations lead- 
ing to a wide variety of different properties. For exam- 
ple, Pr3Ga5SiOi4 has drawn some attention as a two- 
dimensional Kagome system which only magnetically 
orders under chemical pressure J^iiS The only strongly 
magnetic langasite compounds discovered so far are those 
containing iron ions and have been observed to display a 
series of novel piezoelectric and magnetic properties;^ 

A particular iron based langasite compound is 
Ba3NbFe3Si20i4 which consists of Fe^+ S = 5/2 ions. A 
detailed structural study of this compound has been pre- 
sented previously illustrating a unique single phase mag- 
netic spiral structure with each triangle adopting a con- 
ventional 120° arrangement but with a spiral pitch along 
the c axis resulting in a magnetic propagation vector of 
qo=(0, 0,1/7).^^ The underlying structure is displayed in 
Fig. [T] Panel a) illustrates the Fe^^ framework in the 
a — b plane and shows the material consists of triangles 
of magnetic ions arranged on a hexagonal lattice. The 
lower panel b) illustrates the framework along the c axis. 
The space group is chiral and the crystal possess a hand- 
edness, as seen from the "helical" paths along the c-axis 
(panel b). The exchange pathways labeled Ji=i_5 are dis- 
cussed later in this paper. The magnetic order occurs at 
T^r=27 K despite a Curie- Weiss temperature of &w=- 
173 K, owing to the effects of low-dimensionality and 
frustration. Further work on this material has suggested 
the possibility of multiferroic properties at low tempera- 
ture and as a function of applied magnetic field i^°i^^ 

We present a neutron elastic and inelastic scatter- 
ing study of the fluctuations and critical properties in 
Ba3NbFe3Si2 0i4. We show that the magnetic chirality is 
imposed by symmetric Heisenberg exchange only without 
the need for a dominant Dzyaloshinskii-Moriya interac- 
tion and illustrates a strong coupling between structural 
and magnetic chirality. The paper is divided into four 
sections discussing our results along with an introduc- 
tion and a section on the experimental details. We first 
investigate the temperature dependence of the magnetic 
correlations near the critical wave vector through the use 
of inelastic powder and single crystal studies to show 
that the transition behaves classically as expected for a 
two dimensional triangular magnet. This section also in- 
cludes our results using polarized neutrons to investigate 
the energy integrated diffuse scattering as a function of 
temperature. To motivate our inelastic work, we then 
present our results from spin-dimer calculations which 
provide an estimate for the relative values of the mag- 
netic exchange parameters. We then discuss the spin- 
waves and our heuristic model to extract the exchange 
constants. We finish the paper with a derivation of the or- 
dering wave vector and Curie- Weiss temperature. While 
Dzyaloshinskii-Moriya interactions are present owing to 
the symmetry of the lattice, we find that we can un- 
derstand the ordering wave vector (qo=(0,0,~l/7)), the 
spin- waves, and Curie Weiss constant (Qw) in terms 




FIG. 1. (Color online) Crystal structure of Ba3NbFe3Si20i4 
projected in the a6-plane (top panel) and along the c-axis 
(bottom panel). The Fe04 tetrahedra are shown in light blue 
color. The Fe and O atoms are represented as large and small 
spheres respectively. The Ba,Nb and Si atoms are omitted for 
clarity. Exchange interactions are shown in plane by thick line 
(Ji) and thin lines (J2). Along the c-axis, the three inequiva- 
lent exchange interactions are marked by their respective 0-0 
paths, shown as thick bonds. 



of a model based solely on Heisenberg exchange with a 
magnetic chirality introduced because of the handedness 
forced by the crystal structure. 



II. EXPERIMENTAL DETAILS 

Neutron experiments utilized instruments at both the 
ISIS spallation neutron source (Rutherford Appleton 
Laboratory, UK) and the FRM2 reactor (Germany). To 
initially characterize the dynamics and temperature de- 
pendence, we used the MARI direct geometry chopper 
spectrometer at the ISIS facility. The sample consisted 
of 11 g of powder mounted in an annular geometry and 
cooled using a closed cycle displex. A Gd Fermi chop- 
per was used to define an incident energy in parallel with 
a disc chopper to suppress background and high-energy 
neutrons above ~ 250 meV. Fast neutrons, with energies 
in excess of ~ 1 eV, were removed using a nimonic (or to) 
chopper located close to the target face. To measure the 



magnetic fluctuations as a function of energy and tem- 
perature we used a fixed incident energy of Ej=10 meV 
taken with a Fermi chopper speed of 250 Hz. 



The single crystal sample, with mass 6 g, was grown 
using the fioating zone technique and was characterized 
at the test and alignment spectrometer at the ISIS fa- 
cility using a four-circle geometry. The lattice constants 
were measured to be a=b= 8.539 A, and c=5.241 A with 
7=120°. The sample was aligned such that Bragg reflec- 
tions of the form (HOL) lay within the horizontal scat- 
tering plane. Two different sets of experiments were per- 
formed at the FRM2 reactor using a closed cycle refrig- 
erator to reach temperatures as low as 2.5 K. To mea- 
sure the spin-waves across the entire Brillouin zone, the 
PANDA cold triple-axis was used. A vertically focused 
pyrolytic graphite (PG002) monochromator was used to 
select an incident neutron energy refiected onto the sam- 
ple position. A fiat PG002 analyzer was used to select 
the final energy of the neutrons analyzed in the detector. 
A cooled Beryllium filter on the scattered side filtered 
out higher order neutrons refiected off the monochroma- 
tor and coUimation sequence was set to open-SO'-^-SO'- 
open. Final energies of E/=5.0 and 2.5 meV were chosen 
to study the low-energy spin waves and E/=13.5 meV 
was used to study the magnetic order parameter with 
the Beryllium filter replaced by a graphite filter. All 
data taken on PANDA have been corrected for higher 
order contamination of the incident beam monitor which 
determines the counting time. The correction factor is 
discussed in detail elsewhere<2^^ 
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FIG. 2. (Color online) a) The temperature dependence of 
the elastic magnetic scattering measured with a final energy 
of E/=13.5 meV taken on the PANDA cold triple-axis. The 
curves are power laws descrbed in the text. Panels b) — d) are 
scans along the (OOL) direction illustrating that the incom- 
mensurability remains fixed with the temperature. 



III. PARAMAGNETIC SCATTERING AND 
MAGNETIC ORDER 



To measure the energy integrated diffuse magnetic 
scattering, we have used the DNS polarized cold 2 axis 
diffractometer located at FRM2. An incident energy of 
3.64 meV was selected using a horizontally and vertically 
focused PG002 monochromator. The beam was polarized 
using m — 3 Scharpf supermirror polarizers. The polar- 
ization at the sample was fixed through the use of an xyz 
coil with the x direction chosen to be parallel to the av- 
erage Q at the sample position and z vertical. With the 
use of hipping coils in the incident and scattering beams, 
the two spin-fiip and non spin fiip cross sections could be 
measured with the neutron polarization along the three 
orthogonal Cartesian coordinates (a total of 12 channels). 
The flipping ratio was 20 ± 1 and was not found to devi- 
ate from this value regardless of the direction of neutron 
polarization. All spin-flip data has been corrected for 
the feed through from the non-spin-flip channel. The 
scattered beam was measured with 24 detectors equally 
spaced 5° apart covering a total angular range in scatter- 
ing angle (26*) equal to 120°. These measurements were 
performed using a two-axis geometry and therefore pro- 
vided an approximate measure of the energy integrated 
intensity {S{Q)). 



A. Magnetic order and wave vector 

A summary of the elastic scattering which character- 
izes the magnetic structure is presented in Fig. [5] and 
illustrates the temperature dependence of the magnetic 
Bragg peak at (5=(l,0,~l/7). The temperature depen- 
dence of the magnetic Bragg peak is plotted in panel a) 
and was taken using the PANDA cold triple-axis with 
E/=13.5 meV. The solid red curve is a plot to a power 
law I (X {Tc- T)2/3 with Tc=27 K and ^=0.25 follow- 
ing the critical exponents measured in CsMnBrs, and 
CsVCls."^^ All of these systems have demonstrated to be 
XY-like in nature with a 120° structure similar to that 
in one trimer in Ba3NbFe3Si20i4. 

Panels b) — d) in Fig. [2]illustrate scans along the (OOL) 
direction at several temperatures. The results illustrate 
that while there is a significant change in the intensity, 
the wave vector of the magnetic peak does not change 
within error even close to the Neel transition tempera- 
ture. This point is unusual with chiral magnets or in- 
commensurate systems as the chirality (and hence wave 
vector) can be considered as a measure of the order pa- 
rameter and this should vary smoothly near T^v unless 
the transition is first-order. For example, the field driven 
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FIG. 3. (Color online) The three spin-flip cross sections mea- 
sured with the incident beam polarized parallel to x, y, and 
z. A cut along the L direction is displayed in panel d) for all 
three cross sections. The cut was integrated over the range 
0.9 < H < 1.1. 



commensurate-incomniensurate transition in CuGeOa is 
first order and becomes second order on doping with non 
magnetic impurities such as Mg resulting in a contin- 
uous change of the incommensurate wave vector near 
TAr.~ Another example is the incommensurate transi- 
tion in Rb2ZnCl4 which exhibits a strong temperature 
dependence of the wave vector near the critical temper- 
ature.— To track the incommensurability as a function 
of temperature we now investigate the energy integrated 
intensity which is characteristic of the critical spin fluc- 
tuations using polarized neutrons. 



B. Polarized diffuse scattering measurements 

The spin of the neutron allows a sensitive measure of 
both the phase and magnitude of the magnetic cross sec- 
tion. The cross sections for polarized neutrons have been 
studied theoretically and are presented in Refs. l27H30i 
Of particular interest is the spin-flip (SF) cross section 
when the incident beam is polarized parallel to the mo- 
mentum transfer (Q), 



FIG. 4. (Color online) The spin flip cross section with the 
incident polarization parallel to x is illustrated at a) 30 K 
and h) 100 K. The data are plotted on a log scale. The data 
were taken on DNS with an Ei=3.64 meV. 



Here p is the magnetic scattering length and 5^ is the 
value of the magnetic spin perpendicular to the momen- 
tum transfer. The symbols H — / — I- denote the cross 
sections when scattering from a spin-up to a spin-down 
state and from a spin-down to a spin-up respectively. In 
collinear magnets, the values of these cross sections are 
the same as the second term in the above equation is 
identically zero. 

In a chiral magnet, one can define a vector chirality 
via a vector product of two neighboring spins, averaged 
over three spin pairs by 



3\/3 



Et^* X S,], 



(2) 



iz-{S,^y.S,A_)]- (1) 



which is related to the cross product in the cross sec- 
tion written in Eqn. [T] In a chiral magnet the cross 
product in Eqn. [2] is not necessarily equal to zero and 
I"' is not equal to I ^ . This property of chiral magnets 
and polarized neutron scattering was first demonstrated 
in MnSi where the Dyzaloshinski-Moriya interaction re- 
sults in spiral magnetic correlations 4i Given that in our 
notation x is parallel to the average Q, a subtraction 
of scattering channels in -\-x and —x channels gives the 
cross product term above, sensitive to the chirality. The 
cross product in Eqn. [T] does not appear when the inci- 
dent beam polarization is perpendicular to the momen- 
tum transfer (Q) . This is the case of for y and z polar- 
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FIG. 5. (Color online) The spin flip cross section with the 
incident polarization parallel to x and ~x. Panels a) — c) 
illustrate the difference between the spin flip cross sections 
7^ and / ^ at several temperatures. Panel b) and c) are 
linear cut along (I,0,L) illustrating the difference between the 
two cross sections. The cuts was integrated over the range 
0.9<H < f.l. 



izations and these two channels can be used to determine 
the magnetic moment direction or the orientation of S± 
in the above formula. 

Fig. [3] provides an overview of the magnetic scattering 
cross section in the paramagnetic phase by displaying 
the spin-flip intensities for polarizations parallel to the 
three orthogonal directions x, y, and z at T=50 K. A cut 
through the data along the (1,0, L) direction is plotted in 
panel d) for the three polarizations. The spin-flip cross 
section parallel to x is a measure of the entire magnetic 
cross section. As illustrated in panel d), the spin flip cross 
section parallel to x is approximately twice that of the 
y and z channels. This indicates that the scattering is 
isotropic to a first approximation, as expected for para- 
magnetic scattering above T^v- There is a difference be- 
tween the spin-flip cross sections in the y and z channels 
with the z cross section being larger than the y channel at 
Q={lfl,0). This indicates a small anisotropy between the 
paramagnetic scattering above T^r and maybe expected 
given the highly anisotropic structure reported)^ 

Fig. 13] illustrates the temperature dependence of the 
contours of magnetic scattering above Tn at 30 K (panel 
a)) and 100 K (panel b)). Both data sets were taken 
with the incident beam polarization parallel to x and 
therefore are a measure of the total magnetic scatter- 



ing cross section. The data are plotted on a logarithmic 
scale. Since x is parallel to the average Q and no an- 
alyzer crystal was used, these scans provide an approx- 
imate measure of the total energy integrated magnetic 
scattering cross section. At T=30 K, well defined dif- 
fuse scattering contours exist around the incommensu- 
rate positions at <7o=(0,0,~ 0.15) where magnetic order- 
ing was observed below T^r. At T=100 K, clear diffuse 
scattering is present though much reduced in intensity 
and considerably broader in momentum indicating a de- 
crease in the correlation length. It is interesting to note 
that the contours of constant scattering intensity imply 
an anisotropic line shape which is broader in the [100] 
direction than the [001]. This implies that the c-axis 
correlations are significantly stronger and therefore it is 
important to understand the c-axis magnetic exchanges 
in understanding the critical and physical properties of 
this system. We will return to this point later when dis- 
cussing the spin- waves at low-temperatures. 

As pointed out in Ref. [S^, a vector chirality variable 
can be defined by Eqn. [5] and is proportional to the 
cross product in the second term of Eqn. [TJ This chiral- 
ity should occur simultaneously with spin ordering and 
follow a power law divergence similar to the magnetic or- 
der parameters in conventional second order phase tran- 
sitions. This order parameter can be derived in polar- 
ized neutron scattering experiments as the difference be- 
tween the two spin-flip cross sections I^ equal to / ^ 
when the neutron polarization is parallel to the momen- 
tum transfer vector. Therefore, near the phase transition 
the difference between the two spin-flip channels and the 
value of the incommensurability should to go zero. We 
have investigated this quantity using the DNS polarized 
neutron diffractometer. 

The difference between the SF cross sections /"' and 
/ ^ with the incident neutron polarization parallel to 
the x axis is displayed in Fig. [5] Panels a — c) plot 
subtracted intensity /'' — / ^ and is sensitive to the 
chirality parameter discussed above. If there were no 
chirality present, then the cross product in the Eqn. [T] 
would vanish and the difference would be zero. All panels 
at temperatures ranging from 30-100 K clearly show the 
presence of correlated magnetic scattering with a vector 
chiral component, even at 100 K well above the order- 
ing temperature of 27 K. The magnetic scattering is also 
clearly displaced away from the nuclear Bragg peak po- 
sitions and located close to the position in momentum 
where the magnetic Bragg peak is observed at low tem- 
peratures. As these measurements integrate in energy, we 
are not able to determine from these scans if this diffuse 
scattering is dynamic or static or on what time scale the 
fluctuations occur at. This will be discussed in the next 
section with energy resolved triple-axis and direct geom- 
etry chopper measurements. The white gaps near the 
nuclear Bragg peak positions in Fig. [5] are regions where 
the strong feed through from the non spin flip channel as 
a result of nuclear scattering did not subtract well and 
hence have been removed from the data for clarity. 
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FIG. 6. (Color online) The spin flip cross section with the 
incident polarization parallel to z and —z. Panel a) illustrates 
contours of the raw intensity with beam parallel to z and 
hence shows l^~ . Panel b) illustrates a comparison between 
the I^~ and /— +. The cuts was integrated over the range 
0.9 < H < 1.1. Panel c) is a subtraction of the two cross 
sections illustrating no difference within experimental error. 



To confirm the experimental setup, we plot a similar 
analysis but with the incident polarization parallel to the 
z axis in Fig. [6l This is an important check as the dif- 
ference observed above could result from differences in 
the flipping ratios between the +x and —x configura- 
tions. One of the two cross sections (/"' ) is illustrated 
in panel a) with panel b) plotting a cut of both /+~ 
and / ^. Panel c) confirms that there is no difference 
between the two cross sections to within experimental 
error except near the Bragg position (L=0) where the 
feed through from the non spin-flip cross section con- 
taminates the data and does not subtract out fully. This 
complete subtraction is the result expected for the spin- 
flip cross section when the incident beam polarization is 
perpendicular to the momentum transfer. This conflrms 
the quality of the experiment, the results, and the con- 
clusions applied to the data taken with the polarization 
parallel to x described above. 

The results here are quite surprising given that the chi- 
rality (as deflned in Ref. |32| ) is predicted to diverge at 
T^v, yet we still observe a strong chirality at tempera- 
tures more than three times T^r. This is in contrast to 
the results reported in MnSi which is a prototypical spiral 
magnetic resulting from the Dzyaloshinski-Moriya inter- 
action. In that work, it was found that the difference be- 
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FIG. 7. (Color online) A series of contour plots of the mag- 
netic scattering in Ba3NbFe3Si20i4 at a series of tempera- 
tures. The data were taken on the MARI direct geometry 
chopper spectrometer located at ISIS with a fixed incident 
energy E,; = 10 meV and at a Fermi chopper speed of 250 Hz. 



tween the two cross sections became signiflcantly smaller 
at high temperatures above the ordering temperature. 
Indeed a Landau free-energy expansion did show that the 
two cross sections (/^ and /""'") should become equal 
at high temperatures, consistent with the experimental 
data. This is clearly not the case in Ba3NbFe3Si20i4 
and implies that a picture utilizing Dzyaloshinski-Moriya 
must be examined carefully. In the next section we dis- 
cuss the energy resolved critical dynamics measured with 
both a chopper instrument and a cold triple-axis. 



C. Spin fluctuations and critical properties 

To understand the critical properties, we investigated 
the critical spin fluctuations near the ordering tempera- 
ture using both powder and single crystal samples. The 
inelastic spectrum at several temperatures in the powder 
is summarized in Fig. [7]with the data taken on MARI us- 
ing an Ei=10 meV. The low-temperature (T=4 K) spec- 
trum will be discussed in detail later but the higher tem- 
peratures illustrate the presence of strong magnetic flue- 
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FIG. 8. (Color online) A series of constant-Q cuts taken from 
the Ei = 10 meV MARI data above T]v. The solid curves are 
fits to the modified Lorentzian discussed in the text. The 
data illustrate strong critical fluctuations which broaden with 
increasing temperature. 
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tuations above the ordering temperature of TAr=27 K. To 
extract a line width as a function of temperature we have 
fit a series of constant-Q cuts to the following formula, 



S{E)=xoHE) + l] 



ET 



T^ + E^ 



(3) 



where xo is the real part of the susceptibility, T is the 
energy half width, and n{E) + 1 is the Bose factor. This 
formula allows us from a constant-Q scan to extract the 
real part of the susceptibility {xo) and an energy line 
width (r). Examples of the resolution convolved fits ex- 
tracted from the MARI data above T^r are shown in Fig. 

ISl 

Constant-Q cuts taken using the PANDA cold triple- 
axis spectrometer on single crystal samples are presented 
in Fig. [9] These scans were performed at the measured 
ordering wave vector and the solid lines are fits to Eqn. 
12] plus a delta function centered at the E = elastic 
position defined by the energy resolution. Similarly, from 
these scans we are able to extract parameters for the 
susceptibility xo and the energy line width F. 

A summary of the temperature dependence of xo and 
F are illustrated in Fig. [TU] The scaling factor for the 
susceptibility between the MARI and PANDA data were 
obtained from a common temperature taken at 100 K. 
The two data sets overlap very well illustrating a consis- 
tent analysis. A plot of Xo ^ (Fig- UHl »)) illustrates a 



FIG. 9. (Color online) A series of constant-Q scans taken 
on the PANDA cold triple-axis using a fixed final energy of 
E/=5.0 meV at Q=(l, 0,1/7). The solid curves are fits to the 
modified Lorentzian discussed in the text plus a the resolution 
function centered at the elastic energy (E=0 meV). 



smooth function which intercepts the x-axis at the Neel 
temperature indicating that the transition is likely sec- 
ond order. The solid curve is a plot of xo '^ {T — ^at)"^ 
with 7=1.1. This is the same critical exponent predicted 
and measured for the strongly two dimensional triangu- 
lar magnets CsMnBrs, and CsVCls.'^^ We consider these 
as good systems to compare the critical properties of 
Ba3NbFe3Si20i4 to as these triangular magnets display 
120° magnetic ordering and weak coupling between tri- 
angular planes. Also, because of strong dipolar interac- 
tions, these systems also display a strong XY character 
with the spins being forced in a particular plane. 

A plot of the line width 2F as a function of temper- 
ature above T^r is illustrated in panel b) of Fig. [TOl 
The solid line is a plot to the power law F c>c (T — T^Y 
with 1^=0.53 as found for the two dimensional triangular 
systems described above. The line width (F) can be re- 
lated to the correlation length using hydrodynamic the- 
ory which states F oc fic/S^^ where fic is the spin-wave 
velocity, and ^ is the correlation lengthii^ Given the cor- 
relation length scales as ^ ex (T — Tm)~^ we have used 
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FIG. 10. (Color online) a) A plot of the inverse of the real 
part of the susceptibility 1/xo as a function of temperature 
for data taken on both powder (MARI) and single crystal 
(PANDA) samples, b) The full-width (2r) is plotted against 
temperature. The dashed line is 2r = ksT. The solid curves 
in both panels are the power laws described in the text and 
previously measured in two dimensional triangular antiferro- 
magnets. 



the same power law form to describe the line width in 
panel b) of Fig. [TOl 

The dashed line in Fig. [TO] 6) is a plot of 2T = ksT 
and illustrates that the line width we observe is al- 
ways less than the energy scale set by the tempera- 
ture (ksT). This indicates that the dynamics and crit- 
ical fluctuations are in the classical limit in contrast to 
frustrated magnets with lower spin- value, for example 
NiGa2S4 or in the superconductor YBa2Cu06-i-£!;, where 
2r > ksTJ-iS^^ These results indicate that the criti- 
cal scattering in Ba3NbFe3Si20i4 is indeed classical in 
nature. 

In summary, we have combined chopper and triple- axis 
measurements on single crystals and powders to obtain 
the temperature dependence of the energy linewidth (2r) 
and the susceptibility (xo)- The temperature dependence 
is consistent with the power laws derived in other two- 
dimensional triangular magnets where 7 = l.l,/3 = 0.25, 
and v — 0.53 and indicates that the transition is second 
order and that the fluctuations near Tjv behave classi- 
cally^ 



IV. SPIN DIMER CALCULATION 

The crystal structure for the left-handed chirality is 
represented in Fig. [1] projected in the a&-plane (top 
panel) and along the c-axis (bottom panel). Exchange 
paths are indicated by the label Ji (i=l,5) following the 
conventions used in Ref. J2l]. All exchange terms are 
dominated by contributions from super super-exchange 
(SSE) through Fe^+-0-0-Fe^+ since there are no cova- 
lent Fe'^"'"-0-Fe'^"'" bonds. Ji and J2 are respectively the 
nearest-neighbor and next-nearest neighbor exchange in- 
teractions in the a&-plane. Along c, there are three in- 
equivalent exchanges paths connecting adjacent layers. 
J4 connects each of the Fe'^"'' ions in a triangular unit to 
ions directly above, i.e. related by the lattice translation 
(0,0,1). J3 and J5 links different Fe^+ ions of adjacent 
triangles and create helical paths along the c-axis of op- 
posite chirality (anticlockwise for J3 bonds and clockwise 
for J5 bonds). The relative strengths of the different SSE 
terms are determined by the overlap of the O p wave 
functions, which highly depends on the value of the O- 
O distance with respect to the Van der Waals distance 
(2.8OA) and the value of the Fe3+-0-0-Fe3+ dihedral an- 
gle. Relevant bond distances are presented in the upper 
row of Table HVl 

The relative strengths of the five SSE interactions have 
been derived semi-quantitatively by a spin-dimer analysis 
based on extended Huckel tight-binding (EHTB) calcu- 
lations.'^^ The atomic orbitals for Fe and O are approxi- 
mated by double-C Slater orbitals, using the parameters 
listed in the supplementary information of Ref. [SJ- The 
Fe'^+ ion is in a distorted tetrahedral configuration (point 
symmetry 2), and the non-bonding Eg levels and anti- 
bonding T2g levels found for Td symmetry splits into sin- 
glets Eg ->2A and T2g ->A©2B. The molecular orbitals 
of the Fe04~ complex calculated by EHTB are shown in 
Fig. [TTJ Each the the singlet (labeled /x=l,5 from the 
lowest to the highest energy) is occupied with a single 
electron due to the high-spin d^ electronic configuration 
of Fe^+. The antiferromagnetic contribution of each of 
the SSE interaction can be evaluated in turn by consid- 
ering the corresponding Fe20g ~ dimer and calculating 
the average of the squared energy difference between the 
bonding and antibonding states: 



J(x((Ae)^)(x V {Ae,,) 



M=l,5 



(4) 



where the sum runs over the five sets of bonding- 
antibonding states resulting from correlations between 
pairs of all the 3A and 2B molecular orbitals. Essentially 
the values of < (Ae)^ > are proportional to the square 
of the hopping integral (t^), i.e. to the exchange ener- 
gies. The results are listed in Table HVl The strongest ex- 
change interaction is Ji, slightly stronger that the other 
in-plane exchange J2. Along the c-axis, J5 is by far the 
strongest exchange, as already pointed out in Ref. l2ll . 
and approximatively four times larger than J4. Due to 
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FIG. 11. (Color online) Molecular orbitals of the Fe04~ com- 
plex calculated by the extended Huckel tight-binding method 
(see text for details). The five levels (all singlets) are shown 
with their corresponding energy in eV, symmetry label, and 
drawing of the molecular orbitals. On the left, a Fe04 tetra- 
hedra is shown in the same projection along the c-axis. The 
unit-cell axes, and position of the two-fold symmetry axis are 
also shown. 



V. LOW TEMPERATURE SPIN- WAVES 

In this section we discuss the spin-waves measured in 
the Neel ordered state at low temperatures. We first 
present the data measured both along (OOL) and along 
(HOO) directions. We then present a heuristic model for 
the spin-interactions which map onto an XY Hamilto- 
nian consistent with our analysis of the critical fluctua- 
tions which is strongly suggestive of a such a character 
for the spins. This model provides a means of estimating 
the exchange constants along [OOL] (Ji=3_5). We then 
derive estimates for the exchange constants within the 
a — b plane (Ji=i.2). Using the values for the exchange 
constants derived from the spin-waves, we derive expres- 
sions for the the ordering wave vector and the Curie- 
Weiss constant and compare with experiment. We focus 
our discussion on the spin dynamics along the c* direction 
as these have the most relevance for the incommensura- 
bility. 



A. Experimental Data 



the very long 0-0 distance, J3 is negligible in comparison 
with all other exchange terms. 
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Applying Eqn. |4]and normalizing to J4, we summarize 
the following for the values of Ji_5, 



Ji = 6.8 

J2 = 4.2 

J3 = 

J4 = l 

J5 = 3.6 



(5) 



We now investigate the low-temperature spin-waves 
which are characterized by the exchange constants and 
use these values to motivate a spin-wave model to de- 
scribe the dispersion curves. 



In this section we present measurements of the spin- 
waves which characterize the low-temperature magnetic 
ground state. To obtain an overview of the magnetic 
scattering we first conducted preliminary measurements 
on MARL The low temperature spectrum is illustrated 
in Fig. [7] a) which illustrates a total spin band width 
of about 4-5 meV and a mode extending to the lowest 
energies at Q^ 0.85 A~^. This wave vector is consistent 
with value of magnetic wave vector of Qo={l, 0,^1/7) 
derived from magnetic structure. 

To investigate the spin- waves further, we performed 
measurements on a single crystal using the PANDA cold 
triple-axis spectrometer. Given the spiral pitch along 
the (OOL) direction, we focused our measurements on the 
spin-wave dispersion along c* in an effort to determine 
the origin of the incommensurate wave vector. The re- 
sults are illustrated in Fig. [T2]where panel a) illustrates a 
contour plot of the measured intensity and panels 6) — d) 
are representative constant-Q scans. While the spectrum 
is quite complex, it can be understood in terms of three 
branches. There is one branch which originates from the 
incommensurate Bragg peak (Q— (2,0,1/7)) and rises up 
to ~ 2 meV at the nuclear zone centre {Q={200)). At 
Q=:(200) we observe two gapped modes. One disperses 
up to ~ 4.5 meV at the zone boundary {Q={2flfl.5)) and 
the other one disperses only up to ~ 3 meV. The solid 
curves in Fig. [12] b) — d) are fits to Lorentzians without 
resolution convolution to extract the energy position. 

We investigated the low-energy spin-waves near the 
incommensurate ordering wave vector of Q— (2,0,1/7) 
by using a fixed final energy of E/=2.5 meV. As out- 
lined above in the experimental section, this provided in- 
creased energy resolution. A summary of the results are 
displayed in Fig. [T^lwith panel a) illustrating a constant 
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FIG. 12. (Color online) a) Contour plot of the spin-waves 
measured along the (OOL) direction. The data are plotted on 
a logarithmic scale. Panels h) — d) are representative constant- 
Q cuts at a serious of wave vectors throughout the Brillouin 
zone. 



Q scan performed at the ordering wave vector. A con- 
stant energy scan at low energy transfers of E=0.2 meV 
is shown in panel h) . The sohd curves are resolution con- 
volved fits to a gapped mode and a gapless excitation 
originating from the incommensurate Bragg peak. The 
inclusion of resolution in the fits is necessary to achieve 
the high-energy tail of the peak in panel a) . The results 
illustrate the existence of two modes originating from the 
ordering wave vector at Q={2fi,l/1) - one gapless and on 
gapped with an gap value of 0.35 meV. These modes be- 
come degenerate above the gap value within experimental 
resolution. The dispersion curves along the (OOL) direc- 
tion are sensitive to the exchange constants along the c 
direction iUustrated in Fig. [T] and labeled as J3,4,5. A 
model used to derive these parameters and interpret the 
physical nature of these modes is presented in the next 
section. 

Having shown the experimental results for the spin- 
waves along (OOL) direction, we now show the results 
along the (HOO) direction. A summary of a series of 
constant-Q scans taken along the (H, 0,6/7) direction is 
displayed in Fig. [14] with a false contour plot in panel 
a) and constant Q scans in b — d). The scans illustrate 
that the three modes observed at the incommensurate 
order wave vector disperse and become more well sepa- 
rated along the (H, 0,6/7) direction. The dispersion of 
these modes are sensitive to the in plane exchange con- 
stants, namely the exchange within a given triangle (Ji) 
and the exchange between triangles (J2)- Values of these 




0.1 0.2 

(2,0,L) (r.l.u.) 



FIG. 13. a) Constant Q scan at Q=(2,0,0.15) at T=2.5 K and 
E/=2.5 meV. A gapped mode is clearly present at 0.35 meV. 
Panel 6) displays a constant energy scan performed below the 
gapped mode illustrating the presence of magnetic scattering 
below the gapped mode and the presence of a gapless phason 
mode. 



exchange constants will be estimated and are presented 
in the next section. 



B. Spin waves along (OOL)-XY Hamiltonian 

Most spiral magnets have the chirality imposed by the 
Dzyaloshinski-Moriya interaction (characterized by the 
vector D) which appears in the Hamiltonian as a cross 
product of spins {D ■ Si x 82)- While such a model may 
fit the low-energy spin waves described above, we find to 
consistently describe the observed spin-wave band width 
(defined by J) and the incommensurate wave vector (de- 
scribed by the spiral pitch a = arctan(_D/J)) would re- 
quire D ^ J. This is rather unphysical given that D 
is a relativistic correction related to spin-orbit coupling 
and hence expected to be small in comparison to the 
spin-exchange J for a dP high-spin electronic state with 
quenched total orbital momentum (L—0). Only very 
small orbital momentum could be expected due to the 
slightly distorted tetrahedral configuration of Fe'^+04, 
with a slight off-centering of Fe. 

We have taken a different approach to this problem 
motivated by the structure (Fig. [T|), the spin dimer cal- 
culations, and the polarized neutron scattering results 
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FIG. 14. (Color online) a) A false contour plot summarizing 
the constant Q scans along the (l-q,0,6/7) direction taken on 
PANDA with E/=5.0 meV. Representative constant-Q scans 
are illustrated in panels h — d). 



presented earlier. We now investigate whether the physi- 
cal properties, including the incommensurate order wave 
vector, the spin- waves, the Curie- Weiss constant {Qcw), 
and most importantly, the magnetic chirality, can be un- 
derstood in terms of a Heisenberg only framework. The 
Fe'^"'"-0^^ coordination is plotted in Fig. [l]and illustrates 
that all of the exchange paths are super-super exchange 
involving two oxygen ions. The structure along the c-axis 
(Fig. [1] h)) illustrates that the two next nearest neigh- 
bor interactions (J3 and J5) have very different exchange 
paths imposed by the crystal symmetry. We now discuss 
a model for the spin- waves along the c direction. 

To model the spin-waves, we start with the case of 
chains of coupled triangles (Fig. [1]6)) and therefore the 
Hamiltonian for isolated triangles weakly coupled along 
the c-axis. Experimentally, we believe this is a good 
starting point owing to the fact that the spiral prop- 
agation wave vector is along c and also the fact that 
the experimental dispersion within the a ~ b plane are 
substantially less complicated in comparison to the spin- 
wave dispersion along the [001] direction as illustrated in 
comparing Figs. [T2]and[T4l This approach is also corrob- 
orated by the paramagnetic diffuse scattering presented 
earlier which illustrated that the spin correlations along 
[001] were significantly stronger than [100]. Based on 
the coordination paths illustrated in Fig. [T]we pursue a 
model based on one nearest neighbor exchange (J4) and 
two next nearest neighbor (Js^s). We therefore write the 
Heisenberg component of the spin Hamiltonian as. 



Hha — J4:Si ■ S^,l + Js^j • S'. 



H+l 



J50,- • O, 



i+1- 



(6) 



The index i correspond to different triangles along the c- 
axis and the a, &,and c describe the three different spins 
within a given isolated triangle. The interactions J3,4,5 
define the three different exchange pathways between two 
different triangles along a chain discussed above. The 
coupling between a and a spins on different triangles rep- 
resent the nearest neighbor interaction ( J4) and 6, c spins 
correspond to next nearest neighbors (^3,5). 

To obtain a dispersion for linear spin-waves, we have 
imposed a ground state defined by the standard 120° 
structure within a triangle and a spiral along c, defined 
by the pitch angle a. Following previous spin- wave cal- 
culations for triangular antiferromagnets (Ref. 38), the 
transformations are represented by the following for first, 
a spins, 



Qax_ 
Jj - 

ST- 



' cos{ja) — s^ sin(JQ;) 



s^ cos(ja) 



sj sin(ja) 



for b spins. 



S'j'=^ s| cos(ja + 2tt/3) ~ s^ sm{ja + 27r/3) 
S'^y^ s) cos{ja + 2tt/3) + sj sin(ja + 27r/3) 

nbz ^z 



] 3 

and for c spins. 



(7) 



S-f = s| cos(ja - 27r/3) - s^ sin(ja - 27r/3) 
S^y= s) cosija - 27r/3) + s^ sm{ja - 27r/3) 

By symmetrizing the Hamiltonian and writing in terms 
of the new transformed coordinates Si, the Heisenberg 
component of the Hamiltonian takes the following form. 



Hh,. = F («+i + 44+1) + e («+i) (8) 



with 



r= J4 cos(a) + Jz cos(a -I- 27r/3) -t- J5 cos(a — 27r/3) 
Q^h + ,h + ,h. (9) 

This model therefore maps directly onto a Hamiltonian 
with XY symmetry. We note that this approach does 
not directly provide information on the exchange between 
S = 5/2 spins in a triangle which is expected to be the 
strongest interaction. However, the triangle-triangle cou- 
pling described by our model is crucial to understand the 
spin-waves and chirality along the c-axis. 
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(2,0,L) (r.l.u.) 



FIG. 15. (Color online) A plot of the peak position in energy 
as a function of momentum transfer along the c* direction 
(OOL). The solid lines are fits to the spin- wave models dis- 
cussed in the text. The dashed line is a calculation with J3 
fixed to zero and just fitting J4 and J5. The data has been 
symmetrized around L=0 for clarity. 



This approach follows closely the analysis applied to 
Ba2CuGe207i^'^ The main difference in our calculation 
is that we have imposed a chirality to the Hamiltonian 
through two different next nearest exchange constants 
Js^s instead of using a Dzyaloshinskii-Moriya interaction 
as required in Ba2CuGe207. 

Based on this analysis, we can then use the formulae 
presented in Ref. |4l| to obtain the spin-wave dispersion 
along (OOL). We write the following general formula for 
the predicted modes. 



three modes for a given chiral domain. The first is gapless 
and associated with rotations of the spin within the spiral 
plane. The other two modes correspond to transverse 
excitations of the spiral plane and may be gapped due to 
anisotropy.*''!^ 

This model provides a physical understanding of the 
spin- waves and a fit is illustrated in Fig. [13 Fig. [15] dis- 
plays three branches-one which is gapped and centered 
at the nuclear zone center Q— (2,0,0) (colored green) and 
two branches centered around the incommensurate posi- 
tions (colored red and blue). In our formulation above 
in terms of A{L) and B{L) these would correspond to 
branches such that L = L and L = LzLa. The solid cir- 
cles in Fig. [15] are obtained from constant-Q scans and 
the open circles are obtained from constant energy scans 
done with a fixed E/=5.0 meV. This picture illustrates 
that the complicated dispersion measured and plotted in 
Fig. [12] is simplified by the fact that one branch is as- 
sociated with the spin-waves originating from the other 
incommensurate magnetic Bragg peak. 

The branches near the incommensurate points {L — 
L ± a) consist of both a gapped and gapless excitations. 
The gapless mode is the required phason for a spiral 
structure and corresponds to a rigid rotation of the spins 
within the spiral plane. This mode is reproduced from 
our spin- wave analysis with D — 0. The gapped mode lo- 
cated at the incommensurate position can be interpreted 
as the result of fiuctuations of the spiral plane which 
follow the symmetry of the magnetic order. This neces- 
sitates a finite Di. The large gapped mode at the zone 
centre (L = L) can be interpreted as fiuctuations of the 
spiral plane along the c-axis or perpendicular to the spiral 
plane and can be reproduced with a significant value of 
D2 ■ These fiuctuations have the symmetry of the lattice 
and not the magnetic order and therefore appear cen- 
tered around the commensurate positions. Our analysis 
here relies on the existence of two anisotropics associated 
with the gaps at the incommensurate positions and the 
nuclear zone centre. 

A fit to the above model is illustrated in Fig. [TS] (red, 
blue, and green curves) resulting in the following param- 
eters 



E{L)^2SjA{Lr-B{Lr, 



with Aj^ and Bj^ are derived as 



(10) 



A{L)= e + ^^ J^' cos(27rZ) + A 

B{L)= e cos(27rZ) - ^ ~ " -* cos(27rZ), (11) 

with 6 and F defined in Eqn. [3 We have included a 
single-ion anisotropy term Di in a similar manner to that 
used for TbMnOs and LaMnOs to account for gaps in the 
excitation spectrum, particularly at the ordering wave 
vector and at the nuclear zone centrci^:— We expect 



Ji = 0.10 ±0.02 meV 
J5 = 0.33 ± 0.02 meV 
.h = 0.13 ±0.02 meV. 



(12) 



We also extract two single-ion type anisotropics from 
the curves associated with in-plane and out of plane fiuc- 
tuations Di = 0.0040 ±0.0001 meV and D2 = 0.31 ±0.05 
meV respectively. These anisotropics result in a finite 
gap for spin fluctuations. The large anisotropy found for 
out of plane fluctuations confirms the strongly XY nature 
of the magnetic properties as derived from the tempera- 
ture dependence of the critical spin fluctuations and as 
postulated based on the magnetic structure. 
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These parameters are broadly consistent with the spin 
dimer calculations which predict a helical set of exchange 
constants. The magnitudes are generally consistent with 
the calculation except the value of J3 is much larger than 
predicted but is required to provide a good fit to the spin- 
wave dispersion near the Brillouin zone boundary. The 
spin- wave formulation for the excitations along (OOL) are 
quite sensitive to the values of J3, particularly near the 
zone boundary. This is demonstrated in Fig. [15] by the 
dashed line which is a calculation with J3 fixed to be zero 
and J5=0.5 meV. To obtain the relatively flat dispersion 
of the modes near the zone boundary, a relatively large 
value of J3 is required within our formulation. The spin- 
wave calculation is equally sensitive to J4. While we can 
reproduce the fiat nature near the zone boundary and 
the band-width with competing J4 and J5 values only, 
such a model does not fit the initial slope of the modes 
near the magnetic Bragg peaks. This spin-wave model 
therefore provides a sensitive determination of J3,4,5. As 
a check of our analysis, we now use the derived values of 
J3,4,5 to calculate the ordering wave vector and compare 
the value with experiment. 



1. Incommensurate ordering wave vector qo 

The ordering wave vector defined by the pitch a can 
be obtained by calculating the classical minimum of the 
Hamiltonian defined above. The classic exchange energy 
for the s pira l magnetic structure proposed by Marty et 
al. (Ref. uM . independent of any Dzyaloshinskii-Moriya 
interaction, can be written: 



E{a) - 2Z ■^^^i ' "^i+l + "^3'^*' ' "^i+l + ■^^^i ' '^''+1 ^ 



N{Ji cos(a) -I- J5 cos(a + 27r/3) + J3 cos(a - 27r/3)), 

(13) 

where N is the number of cells considered along the c- 
axis, and the sum (i) runs from 1 to N. This is the en- 
ergy for the left-handed crystal structure, and the neg- 
ative triangular chirality. The energy for the positive 
triangular chirality can be obtained by changing 27r/3 
into -27r/3 and vice-versa. The in-plane exchange en- 
ergy is here irrelevant because it does not depend on 
a. The magnetic modulation is considered incommen- 
surate (r ~ i), i.e. the summation runs over an infinite 
number of sites and only normal terms contribute to the 
energy (Umklapp terms sum-up to zero for an incommen- 
surate structure) . The resulting functional of the energy 
is extremely interesting, because it is not invariant by a 
change of sign of a. Changing a to -a is equivalent to 
permuting J3 and J5 and our previous analysis showed 
that these two exchange interaction have very different 
strengths. This demonstrates that for a given triangu- 
lar chirality, a single magnetic chirality is stabilized by 
isotropic exchanges alone, whose strengths are imposed 



by the chiral crystal structure. The stability condition 
for the energy {dE/da—0) leads to the following incom- 
mensurability (defined by tan(a)) in terms of the three 
exchange parameters (J3,4,5) defined above: 



tan(a) = 



\/3(J3 - J5) 
2 J4 - ( J5 + J3) ■ 



(14) 



Substituting the values extracted from the spin-wave dis- 
persion along c* for J3,4,5, we find 



90 



2tt 



= 0.148 ± 0.007 r.Lw. 



(15) 



in very good agreement with the experimental value of 
0.143 r.l.u. For the opposite triangular chirality, the 
value of qo will be reversed. Both structures are com- 
patible with the results of Marty et al. derived from 
neutron polarimetry^. This illustrates a consistent de- 
scription of the incommensurate wave vector and the 
spin-waves in terms of a Heisenberg only model for the 
spin interactions. This result also corroborates the spin- 
wave analysis which implies that all three competing ex- 
change constants are required to describe the dynamics. 
We note that a J4-J3 or a J4-J5 only picture cannot de- 
scribe the incommensurate wave vector. We emphasize 
also that while the Dzyaloshinskii-Moriya interaction is 
present owing to the low crystal symmetry, our analy- 
sis shows that it is not required and would only be a 
small perturbation to the energy asymmetry between J3 
and J5 which determines the ordering wave vector. In 
fact, one can rewrite equation [T2] to include three dif- 
ferent Dzyaloshinskii-Moriya interactions alongside the 
three exchange paths. We consider only the z-component 
of the DM vectors, that couples to the in plane spin struc- 
ture, labeled by DM^ (i=3,4,5). The stability condition 
becomes: 



tan(a) 



\/3( J3 - J5) + 2DM4 + DM3 - DM5 
2J4 - (J5 + ./a) + V^iDAh - DKh) 



.(16) 



It is clear that for Fe^^ , the strengths of the DM terms 
(proportional to —^J) will be at best only of few per- 
cents of the isotropic terms and will not play any role 
in selecting a given magnetic chirality, unless for acci- 
dental degeneracies of the different J terms. The values 
extracted from the spin-wave analysis indicate that this 
is not the case here. 

Since the single magnetic chirality (sign of a) depends 
on the initial choice of the triangular chirality, the ques- 
tion remains about the absolute chirality of the system. 
However, it is important to realize that the magnetic 
structure corresponding to the positive and negative tri- 
angular chiralities are supported by different irreducibles 
representations of the P321 group (totally symmetric ri 
for positive triangular chirality and T2 for negative). The 
presence of single ion anisotropy terms, which are present 
(however small) due to the two-fold symmetry on the 
Fe'^"*' site, will actually select a unique triangular chiral- 
ity, a unique magnetic state will be stabilized. 
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PANDA, T=2.5 K 
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FIG. 16. (Color online) The dispersion curves within the ab 
plane extracted from the analysis described in the text, a) 
illustrates the peak positions along the (H, 0,6/7) direction. 
The dashed line is used to estimate the exchange within a 
triangle (Jo)- The solid lines derived from the general ex- 
pression described in the text for antiferromagnetic exchange. 
b) plots the peak positions around the zone boundary and 
illustrates relatively little dispersion, indicating small next- 
nearest neighbor exchange, c) shows the dispersion along the 
(1-|-H,0,0) direction. The solid curve is a calculation to the 
same formula used in panel a) , but with a finite value of A to 
account for not being exactly at the magnetic zone center. 






S{S+l)5^Jtq^^{hcq) 



(18) 



with S being the nearest neighbor distance. Substitut- 
ing a spin-vifave velocity of 17.5 meV A derived from the 
dashed line in Fig. [16] a) and a bond length d ^ 3.7 
A(Table|IVl), we obtain Ji=1.6± 0.3 meV. 

To derive values for the exchange constant between the 
triangles within the ab plane, we consider couplings (J2 
as indicated in Fig. [ij between triangles and use the 
following expression for an antiferromagnet. 



E{Hf = 43^ (A^ 



t2 • 2 

Jo sm 



{2nH))) . (19) 



Here we have inserted an extra parameter A to take into 
account that we are not exactly at the zone magnetic zone 
center for dispersion curves like those shown in Fig. [16] c) 
or for the presence of anisotropics as discussed above in 
regards to the spin-waves along the (OOL) direction. The 
solid curves in Fig. [16] a) and c) are calculations with the 
above formula with J2=0.45 meV. 

The spin- wave dispersion around the zone boundary is 
illustrated in Fig. [TH] b) and is sensitive to higher order 
(beyond nearest neighbor) interactions between the tri- 
angles. While it is difficult to separate the modes along 
this direction, our results do show that there is relatively 
little dispersion along this direction in comparison to the 
energy values. This indicates a negligible next nearest 
neighbor interaction between the triangles and such a 
scenario is represented by the flat line in Fig. [TO] b) . 

In summary, based on the spin- wave dispersion curves 
in the a — b plane, we extract the following for the ex- 
change interactions 



2. Spin waves within the a — b plane 

To extract information on the exchange parameters 
within the a — b plane, we investigated the spin- wave dis- 
persion along the (11,0,6/7) and (H,0,0) directions and 
along the zone boundary (H,0,0.5). The spin- wave dis- 
persion along (H, 0,6/7) is presented in Fig. [Tn]a). While 
our analysis above for the spin-waves along the (0,0, L) 
direction do not provide direct information on the ex- 
change within a triangle Ji, the initial slope of the low- 
est energy mode along the (11,0,1/7) can be interpreted 
as a first approximation for this exchange parameter. 
Using the second moment sum rule (Ref. |45i ) which 
relates S{Q,E) to the dispersion E{Q) as {E{QY) — 
J dEE^SiQ, E)l J dESiQ, E), 



{E{Qf) = -S{S + 1) ^ J|[l - cos(g • J)], (17) 



and assuming small momentum transfers q close to the 
magnetic Bragg peak and parallel to [100], we write. 



Ji = 1.6 ±0.3 meV 
J2 = 0.35 ± 0.05 raeV. 



(20) 
(21) 



We now compare the exchange parameters obtained to 
the Curie- Weiss temperature previously reported from 
magnetization measurements. 



VI. CURIE- WEISS TEMPERATURE 

Based on our values for the exchange constants, we 
can obtain an estimate for the Curie- Weiss temperature 
{Qcw) using the following formula. 



ks^cw = 



\si.s 



i)E'^- 



(22) 



where the sum is performed over nearest neighbors. The 
above result is derived from a mean-field approach and is 
based on the molecular field produced at a given site by 
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the nearest neighbors and where the Hamiltonian can be 
reduced to a form ^ 5„ • TJ^K i where H^/ is the molecu- 

n 

lar or Weiss field.'*'* Owing to the fact that the magnetic 
structure in Ba3NbFe3Si2 0i4 is quite complex, it is not 
entirely clear how the Hamiltonian can be simplified to 
such a form to provide a direct value for the Curie- Weiss 
constant. Nevertheless, we can estimate the value as fol- 
lows based on the spin- wave analysis described above and 
the number of nearest neighbors for a given Fe'^^ ion, 



2 X ( Ji + J2 + J3 + J4 + J5)/kB 

- 170K. 



(23) 



The terms for J3,4,5 are quite rigorous as there are two 
nearest neighbors. We have chosen to count Ji twice 
and J2 also 2 times as each Fe'^"'" only has two neigh- 
bors coupled by each exchange pathway. The value ob- 
tained above is in reasonable agreement with the mea- 



sured value of & 



'm~-. 



-173 K and -190 K described in 



Refs. Il8|,l20|j and|21|. The comparison confirms the valid- 
ity of our spin- wave analysis and the extracted parame- 
ters. It also confirms our assertion that a helical Heisen- 
berg model is more appropriate to describe the physical 
properties Ba3NbFe3Si20i4 rather than a model invok- 
ing a large Dzyaloshinskii-Moriya interaction term in the 
Hamiltonian. 



VII. CONCLUSIONS 

We have measured the spin fluctuations in 
Ba3NbFe3Si20i4 at both low-temperatures and around 
T^v with the goal of understanding the ordering wave 
vector and the critical properties. We have found the 
critical dynamics are consistent with those of other two 
dimensional triangular magnets with a 120° magnetic 
structure. However, the critical wave vector and chirality 
are independent of temperature and a chiral nature 
of the diffuse scattering remains at temperatures well 
above Tjv . Through an analysis of the structure using a 
spin dimer calculation, we have followed a scheme where 
the exchange constants along the c-axis are helical with 
two different next nearest neighbor couplings. Using 
this and the fact that an XY picture seems appropriate 
we derived the spin-waves for modes propagating along 
the c-axis with the spins confined to the a — h plane and 
have derived the exchange constants. The experimental 
values are reasonable agreement with the spin dimer 
calculation and reproduce the correct magnetic ordering 
wave vector and Curie-Weiss constants within error. 
These results indicate that the structural and magnetic 
chiralities are strongly coupled in Ba3NbFe3Si20i4 and 
that symmetric Heisenberg exchange is required for this. 
These results point to a unique magnetic system where 
the chirality is imposed by the nuclear structure and not 
the Dzyaloshinskii-Moriya interaction. 
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